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Phobos Laser Ranging (PLR) is a concept for a space mission designed to advance tests of rela- 
tivistic gravity in the solar system. PLR's primary objective is to measure the curvature of space 
. . . around the Sun, represented by the Eddington parameter 7, with an accuracy of two parts in 10^, 

I thereby improving today's best result by two orders of magnitude. Other mission goals include 

, measurements of the time-rate-of-change of the gravitational constant, G and of the gravitational 

'n— «' ' inverse square law at 1.5 AU distances-with up to two orders-of-magnitude improvement for each. 

The science parameters will be estimated using laser ranging measurements of the distance between 
an Earth station and an active laser transponder on Phobos capable of reaching mm-level range 
• resolution. A transponder on Phobos sending 0.25 mj, 10 ps pulses at 1 kHz, and receiving asyn- 

^0 ' chronous 1 kHz pulses from earth via a 12 cm aperture will permit links that even at maximum 

range will exceed a photon per second. A total measurement precision of 50 ps demands a few 
hundred photons to average to 1 mm (3.3 ps) range precision. Existing satellite laser ranging (SLR) 
facilities-with appropriate augmentation-may be able to participate in PLR. Since Phobos' orbital 
O ', period is about 8 hours, each observatory is guaranteed visibility of the Phobos instrument every 

<sJ '< Earth day. Given the current technology readiness level, PLR could be started in 2011 for launch 

>' ^ ' in 2016 for 3 years of science operations. We discuss the PLR's science objectives, instrument, and 

{3J[j| mission design. We also present the details of science simulations performed to support the mission's 

I ' . primary objectives. 
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V«D '. I. INTRODUCTION 

0^ ; 

■ Our solar system is a unique laboratory that offers many opportunities to test relativistic gravity. The true renais- 
sance in tests of general relativity began in the 1970s vi^ith major advances in several disciplines-notably microwave 
spacecraft tracking, astrometric observations with very long baseline interferometry (VLBI), and lunar laser ranging 
(LLR). Because of this technological progress our knowledge of gravity improved significantly [1]. 

In particular, analysis of data obtained from radio ranging to the Viking spacecraft determined Eddington's metric 
parameter^ 7 at the level of 1.000 ± 0.002 [4]. Geodetic measurements with VLBI have yielded value of 7 = 0.99983 ± 
0.00045 [5]. LLR, a continuing legacy of the Apollo program, provided improved constraint on the combination of 
5^ ■ parameters 77 = 4/3 — 7 — 3 [6, 7]. The analysis of LLR data [8] constrained this combination as 77 = (4.0 ±4.3) x 10^^, 
leading to an accuracy of ~0.011% in verification of general relativity via precision measurements of the lunar orbit. 
Finally, microwave tracking of the Cassini spacecraft on its approach to Saturn improved the measurement accuracy 
of the parameter 7 to 7 — 1 = (2.1 ± 2.3) x 10~^, thereby reaching the current best accuracy of ^0.002% provided by 
tests of gravity in the solar system [9] (see discussion in [3]). 

While most of the progress in tests of gravity on interplanetary scales was demonstrated using the microwave 
frequencies, laser ranging techniques gained significant momentum [10]. Laser ranging from the Earth to passive 
targets on the lunar surface routinely operates at centimeter level accuracy [7]. Millimeter- level accuracy LLR data 
has been achieved by the Apache-Point Observatory Lunar Laser-ranging Operation (APOLLO) [11, 12]. Over the 
years, LLR has benefited from a number of improvements both in observing technology [10] and data modeling. Today 



o 
o 



X 



^ To describe the accuracy achieved in the solar system experiments, it is useful to refer to the parameterized post-Newtonian (PPN) 
formalism (see discussion in Sec. IV B and [2]). Two parameters are of interest here, the PPN parameters 7 and the values of which 
in general relativity are 7 = /3 = 1. The introduction of 7 and /3 is useful with regard to measurement accuracies [3]. In the PPN 
formalism, the angle of light deflection is proportional to (7 + l)/2, so that astrometric measurements might be used for a precise 
determination of the parameter 7. The parameter /3 contributes to the relativistic perihelion precession of a body's orbit. 
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LLR is a primary technique to perform high-accuracy tests of relativistic gravity including tests of the Equivalence 
Principle (EP), of time- variability in the gravitational constant, and of the inverse-square law of gravity. 

Building on our experience with LLR data analysis, APOLLO's design, construction and operation, and optical 
communications technology development, we investigated the science case, and the instrument and mission design of 
a surface-landed mission to Phobos that would operate a pulsed laser transponder (time-of-flight) capable of achieving 
mm- level accuracies in ranging between Earth and the orbit of Mars. The resulting Phobos Laser Ranging (PLR) 
experiment will build on the success of LLR, but will break the passive lunar paradigm (strong signal attenuation due 
to energy transfer) and extend the effectiveness of this technique to interplanetary scales (Figs. 1, 6). 




FIG. 1: Concept of a laser transponder link between an observatory on Earth and a laser terminal on Phobos. 

At interplanetary distances, active techniques are required to achieve good signal strength (a benefit of energy 
transfer). The development of active laser techniques would extend the accuracies characteristic of passive laser 
ranging to interplanetary distances [13]. Technology is available to conduct such measurements, achieving single- 
photon time resolution measured in tens of picoseconds (ps). One millimeter of range corresponds to 3.3 ps; millimeter 
range precision can be statistically achieved with a few hundred photons in both uplink/downlink directions. For 
comparison, several-meter accuracies have been demonstrated with radio tracking at Mars [14]. 

Interplanetary laser ranging has been recently demonstrated with the MESSENGER (MErcury Surface, Space 
ENvironment, GEochemistry and Ranging)^ spacecraft that successfully established a two-way laser link over 24 
million kilometers of space between the Earth and the spacecraft, achieving 20 cm range accuracy [15, 16]. In addition, 
a successful transmission of hundreds of laser pulses over the distance of nearly 81 million km (0.54 AU) from Earth to 
the Mars Orbiter Laser Altimeter (MOLA)'^, an instrument on the Mars Global Surveyor (MGS) spacecraft in orbit 
about Mars, have demonstrated maturity of laser ranging technologies for interplanetary applications [17]. In fact, the 
MLA and MOLA experiments demonstrated that decimeter interplanetary ranging is within the state of the art and 
can be achieved with modest laser powers and telescope apertures. Achieving mm- level precision over interplanetary 
distances is within reach, thus opening a way to significantly more accurate (several orders of magnitude) tests of 
gravity on solar system scales [8, 18-20]. 

The PLR experiment would conduct investigations in several science areas, including gravitational astrophysics and 
science investigations of Phobos, Mars and the asteroids over a nominal 3 year science mission. Similar analysis of 
achievable science outcomes for various tests of general relativity was conducted for the radio-science experiments with 
the future BepiColombo Mercury planetary orbiter [21, 22] and also suggested for the interplanetary laser ranging 



^ See details on the MESSENGER mission at http://messenger.jhuapl.edu/ 
^ See details on the MOLA instrument at http://mola.gsfc.nasa.gov/ 
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experiments [8, 23-26]. In this paper we focus on highly-accurate laser ranging measurements taken on interplanetary 
scales and their benefits for gravitational physics. 

This paper is organized as follows: Section II describes the science objectives of the PLR experiment and significance 
of anticipated results for gravitational and other fundamental physics. Section III discusses the trade space for 
interplanetary laser ranging. Section IV presents details of our mission trade studies and results of mission simulations. 
Section V presents a technical overview of the PLR experiment and describes the mission and instrument design, 
operations and methodology. We conclude in Section VI with a summary and recommendations. 

II. THE SCIENCE OBJECTIVES OF THE PLR EXPERIMENT 

Recent progress in observational cosmology has challenged Einstein's general theory of relativity as a model for 
gravitation in our universe. From a theoretical standpoint, the challenge is even stronger-if gravity is to be quantized, 
general relativity will have to be modified. The continued inability to merge gravity with quantum mechanics, together 
with the challenges posed by the discovery of dark energy, indicates that the pure tensor gravity of general relativity 
needs modification or aiigmentation. It is believed that new physics beyond general relativity and the Standard Model 
of particles and fields is needed to resolve this issue [3, 27]. 

The kinds of new physics that can solve the problems above typically involve new physical interactions, some of 
which could manifest themselves as violations of the EP, variation of fundamental constants, modification of the 
inverse square law of gravity at various distances, and corrections to the current model of spacetime around massive 
bodies. Each of these manifestations offers an opportunity for experimentation and could lead to a major discovery. 

To respond to these challenges, PLR will push high-precision tests of fundamental gravity into a new regime, by 
using the spacetime in our solar system to explore the physics of the universe. In particular, PLR will improve 
limits on Eddington's parameterized post-Newtonian (PPN) parameter 7 [28], testing the constancy of Newton's 
gravitational constant, G, searching for a new long range interaction via tests of the Yukawa forces on interplanetary 
scales, and conducting tests of the EP. Similarly to its very successful lunar predecessor, PLR could become a major 
solar-system facility to advance research in fundamental physics; it will also benefit the development of interplanetary 
optical communication. 

A. What is the nature of spacetime? Does PPN pcirameter 7 differ from its general relativistic value? 

Modern gravitational research addresses fundamental questions at the intersection of particle physics and 
cosmology-including quantum gravity and the; very ciarly imiverse. This work generated ideas on large extra di- 
mensions, large-distance modification of gravity and brane inflation in string theory, all leading to experimentally- 
testable explanations for the quantum stability of the weak interaction scale. Recent extensions of gravitational 
models including brane- world models and efforts to modify models of gravity on large scales motivate new searches 
for experimental signatures of small deviations from general relativity on a variety of scales, including solar system 
distances (see, e.g., [3]). 

Given the immense challenge posed by the unexpected discovery of the accelerated expansion of the imiverse, it is 
important to explore every opportunity to explain and probe the underlying physics. Theoretical efforts in this area 
offer a rich spectrum of new ideas that can be tested by experiment on solar system scales [27]. 

The Eddington parameter, 7, whose value in general relativity is unity, is the most fundamental parameter in that 
i(l — 7) is a measure, for example, of the fractional strength of the scalar gravity interaction in scalar-tensor theories 
of gravity [29-31]. Specifically, the quantity 5(1—7) defines corrections to the spacetime around massive bodies. To 
date, the most precise value for this parameter, 7 — 1 = (2.1 ± 2.3) x 10~^, was obtained using microwave tracking to 
the Cassini spacecraft [9] during a solar conjunction experiment. This accuracy approaches the region where multiple 
tensor-scalar gravity models, consistent with recent cosmological observations [32], predict a lower bound for the 
present value of this parameter at the level of (1 — 7) ~ 10~^ — 10~^. Therefore, improving the measurement of 7 
would provide crucial information to separate modern scalar-tensor theories of gravity from general relativity, probe 
possible ways for gravity quantization, and test modern theories of cosmological evolution [27]. (In addition, any 
experiment pushing the present upper bound on another Eddington parameter /3, i.e. /3 — 1 = (1.0 ± 1.1) x 10~^ from 
[7, 8], will also be of interest.) 

Figure 2 shows the anticipated scientific discovery space for PLR in measuring the Eddington parameter 7, as 
predicted by a number of modern theories of gravity and gravity augmentation models [27], presented here as an 
example. In particular, prediction #1 comes from the string-inspired scalar-tensor theories that were developed in 
[30, 33], the limit #2 is from the multi-scalar model [34, 35], #3 is from recent work based on the analysis of scalar- 
tensor gravities done by [31], #4 is from the theory of large extra dimensions [36], #5 is from brane- world theories 
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FIG. 2: Scientific discovery space with respect to the PPN parameter 7 and search for time-variations of the gravitational 
constant as predicted by several modern theories of gravity (indicated here by numbers #1-8, see details in the text and in 
Refs. fl, 271). 



and IR-modified gravity [37, 38], #6 are from the /(i?) gravity models [39], # 7 is from relativistic generalization of 
Milgrom's Modified Neviftonian Dynamics (MOND) theory as developed in [40], and, and finally, the limit # 8 is from 
the scalar-tensor- vector theory proposed by [41]. 



TABLE I: Estimated uncertainties are given in the last three columns for parameters of interest as a function of Phobos lander 
mission duration, with 1 mm laser ranging once per day with 2° SEP cut-off and 67 asteroid mass parameters estimated. Note 
that data for the gravitational inverse- square law are given for the shortest Earth-Mars ranges of 1.5 AU. (The current best 
value of solar J® is from [42], which reported J® = (1.82 ±0.47) x 10~^ . Note that data from the Venus Express mission were 
downweighted in our simulations, which resulted in less accurate estimates for J2.) 



Relativistic Effect 
to be studied by PLR 


Current best 
uncertainty 


Mission duration / N of conjunctions 


1.5 yr / 1 cnj 


3 yr / 2 cnj 


6 yr / 3 cnj 


The Eddington parameter 7 
The Eddington parameter /3 
Strong Equivalence Principle, 77 
Solar oblateness, J2 
Parameter G/G, yr~^ 
Gravitational inverse-square law 


2.3 X 10"" 
1.1 X 10"* 
4.3 X lO""* 
4.7 X 10~® 
7 X 10"" 
2 X 10"® 


3.1 X 10"' 
4.3 X lO""* 
1.5 X 10"^ 
6.9 X 10"* 
1.7 X 10"" 
4 X 10"" 


1.4 X 10"' 
1.6 X 10"" 
2.8 X 10"" 
3.2 X 10"** 
2.8 X 10"^" 
2 X 10"" 


7.9 X 10"** 
9.4 X lO"'^ 
8.8 X 10"^ 
2.3 X 10"* 
1.0 X 10"" 
1 X 10"" 



Table I shows the results of parameter estimation arising from simulated laser ranging to Phobos over 1-6 years of 
operation based on daily 1 mm accuracy range points (see details in Sec. IV D 2). Estimated parameters include orbital 
elements and masses of the planets and Phobos, up to 67 individual asteroids, and densities for 3 classes incorporating 
another 230 asteroids. Based on 3 yr simulations, PLR will be able to measure the PPN parameters 7 and /3 with 
accuracy levels of 2 x 10"^ and 2 x 10"*, respectively, to measure the solar oblateness coefficient J2 with an accuracy 
of 3 X 10"®, to measure the fractional time-rate-of-change of the gravitational constant to 3 x 10"^^ yr~^, and to test 
the gravitational inverse square law to an accuracy of 2 x 10"^^ at ~1.5 AU distances and the strong EP parameter to 
3 X 10~*. Thus, PLR could directly verify important predictions of modern theories with an unprecedented accuracy, 
enabling strong advances in the tests of general relativity. More simulation details will be discussed in Sec. IV D 2. 



B. Do the fundamental constants of Nature vary with space and time? 



The possibility that fundamental physical parameters may vary with space and time has been revisited with the 
advent of models unifying the forces of nature based on the symmetry properties of extra dimensions, such as Kaluza- 
Klein-inspired theories, Brans-Dicke theory, and supersymmetry models. Modern gravitational models that attempt 
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to "complete" general relativity at very short distances (embedding it in a more powerful theory capable of addressing 
phenomena on that scale) or "extend" it on very large distances {r-^ 10^^ cm) typically include cosmologically evolving 
scalar fields that lead to variability of the fundamental constants. 

A variation of the cosmological scale factor with epoch could lead to temporal or spatial variation of the gravitational 
constant, G. A possible variation of G could be related to the expansion of the universe depending on the cosmological 
model considered. Variability in G can be tested in space with a much greater precision than on Earth [6]. For example, 
a decreasing gravitational constant, G, coupled with angular momentum conservation is expected to increase a planet's 
semimajor axis, a, as a/a = —G/G. The corresponding change in orbital longitude grows quadratically with time, 
providing for strong sensitivity to the effect of G/G. 

Lunar and planetary ranging measurements currently lead the search for very small spatial or temporal gradients 
in the value of G [6, 43]. Recent analysis of LLR data [6, 8] strongly limits such variations and constrains a local 
(^1 AU) scale expansion of the solar system as d/a = —G/G — —(6 ± 7) x 10^^'^ yi'^^, including those stemming 
from cosmological effects. Interestingly, the achieved accuracy in G/G implies that, if this rate is representative of 
our cosmic history, then G has changed by <1% over the 13.7 Gyr age of the universe. The ever-expanding LLR data 
set and its increasing accuracy will lead to further improvements in the search for spatial or temporal variability of 
G. PLR represents the next advance in testing for variations in G, see Figure 2. 

With 1-mm ranging accuracy, PLR data will permit a measurement of G/G to ^ 3 x 10^^^ yr^^ in about three 
years, see Table L This anticipated accuracy is a result of our simulations that account for the range perturbations 
due to 67 large asteroids. 

C. Do extra dimensions or other new physics alter the inverse square law (ISL)? 

Many modern theories of gravity, including string, supersymmetry, and brane-world theories, suggest that new 
physical interactions will appear at various ranges thereby modifying Newton's gravitational ISL. While most attention 
has focused on the behavior of gravity at short distances (for a review of experiments, see [44]), it is possible that 
non-compact extra dimensions could produce small deviations from the ISL on interplanetary scales [27, 45, 46]. An 
experimental confirmation of new fundamental forces would provide an important insight into the physics beyond the 
Standard Model. 

tcr* 
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FIG. 3: Limits on the ISL violations [44]; PLR will improve ISL tests in outer solar system. Parameters a and A are respectively 
the strength and range of a new physical interaction. 
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In its simplest versions, a new interaction (or a fifth force) would arise from the exchange of a light boson coupled to 
matter with a strength comparable to gravity. Planck-scale physics could give origin to such an interaction in a variety 
of ways, thus yielding a Yukawa-type modification in the interaction energy between point-like masses. In general, 
the interaction energy, V{r), between two point masses mi and m2 can be expressed in terms of the gravitational 
interaction as 

^\''') = [1 + ae '^'^j, (1) 

where r = |r2 — ri| is the distance between the masses, Goo is the gravitational coupling for r — >■ oo and a and A are 
respectively the strength and range of the new physical interaction. Naturally, Goo has to be identified with Newton's 
gravitational constant and the gravitational coupling becomes dependent on r [1]. Indeed, the force associated with 
Eq. (1) is given by: F(r) = VV(r) = — G(r)mim2 r/r^, where effective gravitational constant G{r) has the form 
G(r) =Goo [1 + « (1 + r/A)e-'-A] . 

By far the most stringent constraints on violation of the ISL at large distances to date come from very precise 
measurements of the lunar orbit about the Earth. Analysis of the LLR data tests the gravitational ISL to 3 x 10~^^ 
of the gravitational field strength on scales of the 385,000 km Earth-moon distance [47]. New LLR facilities [11] 
will reach sensitivity of 10~^^ at Earth-moon distances, while mm-level laser ranging to Phobos could improve the 
corresponding ISL solar system tests by two orders of magnitude at interplanetary scales (see Fig. 3). 



D. Is the Equivalence Principle exact? 

Einstein's Equivalence Principle (EP) lies at the foundation of the general theory of relativity; testing this funda- 
mental assumption with the highest possible sensitivity is clearly important, particularly since it is expected that the 
EP will not hold in quantum theories of gravity. The Equivalence Principle can be split into two parts: the weak 
equivalence principle (WEP) tests the sensitivity to composition and the strong equivalence principle (SEP) checks 
the dependence on mass. In its strong form (the SEP) the EP covers the gravitational properties resiilting from 
gravitational energy itself, thus involving an assumption about the non- linear property of gravitation [7] . In the SEP 
case, the relevant test body differences are the fractional contributions to their masses by gravitational self-energy. 
Because of the extreme weakness of gravity, a test of the SEP requires bodies with astronomical sizes. Although gen- 
eral relativity assumes that the SEP is exact, many modern theories of gravity typically violate the SEP by including 
new fields of matter, notably scalar fields [29, 30, 33, 48]. 

Currently, the Earth-moon-Sun system provides the best solar system arena for testing the SEP. LLR experiments 
involve bouncing laser beams off retro-reflector arrays placed on the moon during the period 1969 to 1973 [7]. A 
violation of the EP leads to a radial perturbation of the lunar orbit that, for the SEP, is Sr (4/3 — 3 — 7) 13 cos D 
meters [49-52], with (3 and 7 being two PPN parameters. LLR investigates the SEP by looking for a displacement 
of the lunar orbit along the direction to the Sun. There are laboratory investigations of the WEP (at University of 
Washington) which are about as accurate as LLR [53, 54]. Recent solutions using LLR data yield a test of the EP of 
A{Mg/Mi)ep = (-0.8 ± 1.3) X 10"^^ (corrected for solar radiation pressure as suggested in [55]), where A{Mg/Mi) 
signifies the difference between gravitational-to-inertial mass ratios for the Earth and the moon. After adjusting for 
laboratory WEP tests [52-54], this result corresponds to a test of the SEP of A{Ma/Mi)sEP = (-1.8 ± 1.9) x lO'^^ 
with the SEP violation parameter 77 = 4/3 — 7 — 3 found to be = (4.0 ±4.3) x 10^'' [6 8] . Further improvements may 
be achieved via mm-precision ranges to the moon [8, 11] and new laser ranging instruments deployed on the moon or 
other planetary surfaces. 

A lander on Phobos can push these tests a bit further. In particular, the Sun-Earth-Phobos-Jupiter system tests 
the SEP in a qualitatively different way from LLR [56]. The SEP polarization effect is much larger for the Earth and 
Mars orbits than for the lunar orbit, but the optimum time span is longer than the nominal 3 yr mission. With 1 mm 
precision ranging, the EP-violating polarization toward Jupiter allows the SEP violating parameter r] to be measured 
to 3 x 10~^ for observations spanning up to three years. Additional improvement will be possible with an extended 
duration PLR mission. 

In summary, the proposed PLR mission will provide tests of relativistic gravity in the solar system to an unprece- 
dented precision. It will test the weak-gravity-and-small-speed regime of the cosmologically motivated theories that 
explain the small acceleration rate of the Universe (a.k.a. dark energy) via modiflcation of gravity at cosmological 
scales. PLR would search for a cosmologically-cvolved scalar field that is predicted by modern theories of quantum 
gravity and cosmology, and also by superstring and brane- world models [3, 27]. With an accuracy of two parts in 
ten million anticipated from PLR measuring the Eddington parameter 7, this mission could discover a violation or 
extension of general relativity, and/or reveal the presence of any additional long range interaction. 
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III. TRADE SPACE FOR INTERPLANETARY LASER RANGING 



Highly-accurate time-series of the travel times of laser pulses between an observatory on the Earth and an optical 
transponder on a solar system object could lead to major advances in several science areas, including tests of gen- 
eral relativity, solar system dynamics, and studies of the target body. Any solar-system object with a transparent 
atmosphere and a solid surface would be a suitable platform to deploy a laser ranging instrument [23] . 

The main objective of our study was to identify a solar system location that would enable a highly accurate 
measurement of the gravitational time delay (given by Eq. (5)), which is our primary science signal. Such a location 
should provide the conditions necessary for a much improved determination of the PPN parameter 7 together with a 
number of other relativistic gravity parameters, as shown in Table I. 

The planet Mars is one of those solar-system bodies that are well suited for gravitational physics experiments. 
Proximity to the Earth and the maximization of signal strength is certainly a major factor in choosing an Earth-Mars 
link. Furthermore, the recent increase of space traffic on the Earth-Mars route makes this location quite feasible for 
this purpose. Therefore, it is quite natural to consider Mars as a location of the remote laser terminal as part of an 
interplanetary laser ranging experiment that has the potential to improve the values of several relativistic parameters 
[8, 56]. However, there are factors that make this location less desirable compared to others. 

As part of the effort to optimize the mission concept, a trade space of laser ranging locations was explored. The 
search was limited to celestial bodies within the solar system, since single element solar orbiting "free-flyer" spacecraft 
were considered outside the desired mission cost category due to the perceived need for their drag-free characteristics 
with attendant high expense. The trade space was bounded by the basic concept of laser ranging between Earth and 
another planetary location that could offer opportunities for solar conjunctions with laser beams passing near the sun 
giving good 7 sensitivity. Mission concepts outside this fundamental construct were not included in the analysis. 

Guided by our initial expectations, we started our studies with Mars being the target body for the experiment. 
However, in recent studies we investigated a trade space of mission alternatives, as opposed to an effort focusing 
exclusively on a narrow refinement of the original mission concept. This approach bore two major results. First, one 
of the alternatives explored proved superior to the original concept. Since the Phobos-based mission proved to be 
equivalent in science and lower in both cost and risk compared to the original Mars-based mission, it became the 
preferred solution. The second result of the trade space exploration was a strong understanding of the suitability of 
major alternative locations within the solar system for ranging-based fundamental gravity science. This provided a 
solid, defensible basis for the conclusion that Phobos is the optimal planetary location for laser ranging to and from 
Earth, given the targeted scientific goals and mission class. 
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FIG. 4: Interplanetary laser ranging trade space explored in the study. 



The high-level results of the trade space exploration are shown in Figure 4. The trade space included Mercury 
(both orbiters and landers). Mars (orbiters and landers) and asteroids. Although Mercury offered a competitive 
science return, the practical challenges of surviving in its environment for any meaningful duration made it infeasible, 
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or at least too expensive to remain viable; also, time and cost of orbit transfer were more significant compared to 
other location alternatives. Due to significant influence of non-gravitational forces affecting the spacecraft motion 
(solar radiation pressure, thrustcr firing, etc.) orbital candidates suffered orbit precision difficulties. Suitable asteroid 
candidates could not be found that exhibited enough science-rich near-sun conjunctions over reasonable mission 
lifetimes to stay competitive. The two winners warranting more investigation were the original Mars surface location 
and Phobos. Phobos was preferred to Deimos as a site for laser ranging: because of its faster orbit, shorter precession 
period (Phobos' orbital precession period is 2.26 yr, while for Deimos this period is 54 yr), and smaller area-to-mass 
ratio. 

Below we present results of a further down-select that resulted in choosing Phobos as the preffered location for the 
laser ranging mission. 

A. Tradeoffs between Phobos and Mars 

Both Mars and Phobos laser ranging missions were studied in detail; overall Phobos is a better choice. Mars 
is dynamically well suited to studying gravitational physics, but ranging through its dusty atmosphere presents a 
challenge and solar panels must be large enough to compensate for dust accumulation. The 8 hours (hr) Phobos orbit 
and rotation periods offer three tracking opportunities per day. Our simulations show (see Sec. IV) that ranging to 
Phobos has good sensitivity to improve the values of several gravitational physics parameters. Phobos edged out Mars 
as the preferred solution due to a supcirior science/cost ratio at lower mission risk. The extra parametcirs needed to 
model Phobos' orbit and orientation accurately are no more daunting than dealing with Mars orientation. Detailed 
comparisons between the two concepts are provided as follows: 

• Phobos landing is not complicated by atmospheric entry and descent, so Entry, Descent, and Landing (EDL) is 
simpler, but landing consequently needs more AV. It is approximately 2 km/sec more using multiple burns to 
match orbits. Mars landings are familiar but risky. 

• With current technology. Mars landing sites must be low altitude since atmospheric breaking is used. That 
shifts latitudes north causing greater seasonal variations of sunlight and temperature, plus longer nights and 
lower Earth elevation angles in winter. Phobos landing sites can be equatorial with less seasonal variation and 
better seasonal elevation angle variation of Earth and Sun. 

• Mars has windblown dust that scatters light making the sky a bright noise source while dimming the signal, and 
dust settles on optics and solar panels. There is no dusty atmosphere on Phobos, but electrostatically levitated 
dust is possible, though less trouble. The Phobos sky would be dark with no absorption. 

• A spacecraft on Phobos would likely last longer than on Mars, since it would not experience dust on the solar 
panels, cold winter nights, and reduced winter power from lower elevation of the winter Sun. 

• Ranging to a lander on Mars is sensitive to precession, nutation, UTf* and polar motion, which orient Mars 
[57-59]. Precession is sensitive to moment of inertia, and nutation is sensitive to interior structure. UTl and 
polar motion are affected by the changing polar caps. UTl is also sensitive to tides. 

• Phobos orbit and physical librations add dynamical complexity to the range model, but those complexities 

trade off with Mars UTl, polar motion, nutations, and geocenter motion. Mars UT/PM has been dealt with 
during previous lander mission, while Phobos physical libration has barely been detected and lacks detailed 
measurement. 

• The Phobos orbit has a strong tidal acceleration in longitude, but this should be tractable. Phobos would be 

affected by Mars gravity, but the gravity field has been determined by orbiting spacecraft. It is benficial to orbit 
Phobos on arrival in order to map its gravity field. Phobos has a much larger area/mass ratio than Mars which 
would make it much more vulnerable to impact impulses. 

• Since Phobos makes three orbits per day, a single station on the Earth can always track during one, and 

sometimes two, orbits. By contrast, the similarity of Mars and Earth rotation periods interferes with daily 
tracking from one Earth site requiring at least two tracking sites at different Earth longitudes for daily coverage. 



UTl is the principal form of Universal Time (UT), see details on various timescales at http://en.wikipedia.org/wiki/Universal_Time 
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• 4 hr nights on Phobos, 1/3 of Mars, may require less stored power. But there are three times as many temperature 

cycles. 

• For radio relay through Mars orbiting spacecraft, a Phobos lander would need to be on the hemisphere facing 
Mars. During part of the Mars year Phobos experiences occupations by Mars as seen from Earth and the Sun 
experiences eclipses by Mars. These occultations and eclipses can be up to 54 minutes long so that as little 
as 2.9 hours per orbit, about 9 hours/Earth day, may be available for ranging and charging batteries, and the 
missing times may be for Earth and Sun high in the sky depending on landing longitude. There would be good 
coverage when the Earth and Sun are furthest north or south such as local winter; the eclipses and occultations 
occur for |Mars referenced declinations| < 21° (Mars obliquity is 25°). There would be increased thermal cycling 
during the two eclipse seasons each Mars year. Eclipses and occultations affect only the Mars facing hemisphere 
of Phobos. 

• There is no Phobos atmospheric range correction required and thus a meteorological package is unnecessary. 

• It is easier to support a line of sight close to the sun from the surface of Phobos (due to lack of atmospheric 
scintillation) than from Mars-a clear advantage in accessing the region from which sensitivity to the primary 
PPN 7 signal is greatest. 

• Communication via radio relay is likely less efficient on Phobos than on Mars due to relay directionality, so 
direct communication with Earth may be preferred. Both sides of Phobos (i.e., Mars- facing and far-side) could 
be used for this purpose. 

• Phobos was preferred to Deimos as a site for laser ranging because of its faster orbit, shorter precession period, 
and smaller area-to-mass ratio. 

The considerations above resulted in our choice of Phobos as opposed to Mars as our primary target body for the 
interplanetary laser ranging experiment. 

IV. PHOBOS LASER RANGING MISSION SIMULATIONS 

In this section we discuss the details of the simulations performed to support the PLR's primary science objectives 
presented in Sec. II. 

A. Design of the sensitivity studies 

In order to assess the potential science return from laser ranging to a suitably equipped lander, numerical simulations 
were performed. The simulations were based on the Solar System Data Processing Software developed at JPL and 
used for the development of planetary cphcmerides for NASA missions. The planetary ephemeris (the most recent 
of which is DE421 [60]) is based on a large variety of measurements including radio range measurements to Mars 
landers (Viking, Pathfinder) and orbiters (MGS, Odyssey, MRO). The software includes the capability to model and 
accurately estimate the measured round-trip distance and its dependence on a wide variety of dynamical parameters. 

The orbits of the planets and asteroids are integrated based on the Einstein-Infeld-HofFmann equations of motion 
derived for the PPN metric for general relativity [61, 62] (see Eq. 2). The orbits are integrated based on nominal values 
for the orbital elements of the planets. Sun, moon, and asteroids, their mass parameters, and PPN parameters. The 
sensitivities to these parameters are evaluated by re- integrating the equations of motion many times with perturbations 
of the initial conditions by small amounts from the nominal for each parameter taken one at a time. This procedure 
generates partial derivatives from small finite differences of orbits with respect to initial conditions, masses, PPN 
parameters, etc. Estimates of the dynamical parameters are made using a wide variety of observational data in a 
global least-squares fit (e.g. [60]). 

B. The PPN equations of motion 

Generalizing on a phenomenological parameterization of the gravitational metric tensor field, which Eddington 
originally developed for a special case, a method called the PPN formalism has been developed (see [2, 63] for details). 

This method represents the gravity tensor's potentials for slowly moving bodies and weak inter-body gravity and is 
valid for a broad class of metric theories, including general relativity as a unique case. The several parameters in the 
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PPN metric expansion vary from theory to theory, and they are individually associated with various symmetries and 
invariance properties of the underlying theory (see [2] for details) . 

If (for the sake of simphcity) one assumes that Lorentz invariance, local position invariance and total momen- 
tum conservation hold, the point-mass Newtonian and relativistic perturbative accelerations in the solar system's 
barycentric frame has the following form: 
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where the indices j and fc refer to the N bodies and where fc includes body i, whose motion is being investigated, [ij is 
the gravitational constant for body j given as \Xj = Gnij, where G is the universal Newtonian gravitational constant 
and jUj is the isolated rest mass of a body j. In addition, the vector is the solar system barycentric radius- vector of 



this body, the vector Vij = rj — is the vector directed from body i to body j, and rj 



■ r,; I . The dimensionless 



parameters 7 and (3 are the two Eddington PPN parameters. General relativity, when analyzed in standard PPN 
gauge, gives 7 = /? = 1; other theories yield different values of these parameters [2]. 

To account for a possible violation of the EP and temporal changes in the gravitational constant one can add 
another perturbative acceleration term ^i^^: 



E 



(3) 



where parameter {[mG/mj]i — 1) signifies a possible inequality between the gravitational and inertial masses that is 
needed to facilitate investigation of a possible violation of the EP (see Sec. 11 D). In addition, Eq. (3) also includes pa- 
rameter G/G, which is needed to investigate a possible temporal variation in the gravitational constant (see Sec. 11 B). 
Note that in general relativity 5f j = 0. 

The accuracy of the relativity tests with interplanetary laser ranging will depend on our knowledge of the solar 
gravity field. A source of uncertainty in these tests is the dimensionless solar quadrupole moment J20. The value 
of J20 from solar oscillation data has been reported as (1.7 ± 0.4) x 10~^ by [64] and as about 10% less than this 
value by [65] derived from the oblateness of the photosphere. These results are consistent with the value expected for 
nearly uniform rotation of the Sun. However, solar oscillation data suggest that most of the Sun rotates slightly slower 
than the surface except possibly for a more rapidly rotating core [66, 67] suggesting the value of (2.0 ± 0.4) x 10"'^. 
Recently Fienga et al. reported determination of J2 at the level of (1.82 ± 0.47) x 10~^ using interplanetary ranges 
including Venus Express [42]. Thus the solar oblateness is not an issue for the laser ranging tests of gravity with a 
lander on Phobos/Mars (although J® may be important for the case of ranging to a Mercury lander). The accuracy 
achievable from interplanetary laser ranging may improve upon these results. For this reason, and because of the 
expected correlation of J2 with ^ and other parameters, it is necessary to include J2 in the solution parameter list. 
The acceleration of a test body in the field of the sun must then be modified to read 
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where Rq is the radius of the Sun and is the latitude of the planet with respect to the solar equator. The terms rj 
and Svj are given by Eqs. (2) and (3) correspondingly. 

To determine the orbits of the planets and spacecraft one must also describe propagation of electro-magnetic 
signals between the two points in space. The corresponding one-way light-time equation, with accuracy sufiicient for 
the purposes of this paper, can be presented as below 
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where ti refers to the signal transmission time at body 1, and t2 refers to the reception time at body 2. r\ = \r\{ti)\ 
and r2 = |r2(^2)| are the positions of the transmitter and receiver with respect to the i-th body (including the sun) at 
transmission and reception times correspondingly, and ri2 = |r2(f2) is their spatial separation. The logarithm 

term in Eq. (5) provides our principal sensitivity to parameter 7. We emphasize that the measurement of this PPN 
parameter constitutes the primary objective of the proposed mission (see Sec. II A). 

The PPN expansion serves as a useful framework to test relativistic gravitation in the context of the gravitational 
experiments. The main properties of the PPN metric tensor are well established and are widely used in modern 
astronomical practice [2, 68 70]. The equations of motion Eq. (4) are used to produce numerical codes for the 
purposes of constructing the solar system's ephcmcrides, spacecraft orbit determination [62, 70], and analysis of the 
gravitational experiments in the solar system [2, 3]. In our analysis we also used this set of equations of motion. 

C. Estimated parameters 

In this section, we discuss a series of studies designed to illustrate the results that might be obtained with an 
interplanetary laser ranging (ILR) transponder on Mercury, Mars and Phobos. We used the same data set as in 
DE421 [60], except for the LLR data and planetary transits. The set of estimated parameters included the set of 
seven parameters, namely the value of the astronomical unit (AU), the PPN parameters 7 and /3, solar J2Q, the time- 
rate-of-change in the gravitational constant G/G, GMEarth/GMmoonj and the EP parameter ([rriG/m/jo — 1) ~ tjUq, 
with Uq = —3.52 X 10~^ being the gravitational binding energy of the Sun [56, 71]. For inverse-square law simulations, 
a Yukawa amplitude was estimated for a fixed length. 

TABLE II: Typical set of solar-system parameters estimated in the study (in addition to those listed in Table I). Note that in 
addition to 8 planetary bodies, the orbital elements for I'lutcj were also estimated. 



Estimated parameters 


Number 


Planetary orbital elements 


54 


Astronomical unit (AU) 


1 


GAfEarth/ GMmoon 


1 


Asteroid GMs 


67 


Asteroid class densities 


3 


Spacecraft biases 


8 


Solar corona corrections 


4 


Mercury shape 


8 


Venus radius 


1 


Mars Lander locations 


9 



Total number: 156 



Our studies combined the current solar-system data set with the simulated ILR data in a simultaneous covariance 
analysis. In addition to the six Table I parameters of interest to gravitational physics. Table II shows a summary of 
the parameters in our model of the solar system. The total set of the estimated solar-system parameters included 
162 parameters. The 8 spacecraft biases used in the study and reported in Table II are for the two Viking landers 
1 and 2 (VLl, VL2), Mars Pathfinder (MPF), Mars Global Surveyor (MGS), Odyssey (ODY), Mars Reconnaissance 
Orbiter (MRO), Mars Express (MEX), and Venus Express (VEX). In addition, 4 parameters used to estimate the 
solar corona corrections came from the analysis of the radio-metric tracking data received from MGS, ODY, MRO, 
MEX spacecraft. Note that data from Venus Express mission were downweighted in our simulations, which resulted 
in less accurate estimates for J2 presented in Tables I, III- VI. 

Initial simulations were done for a laser transponder on a Mars lander. The range from an Earth station to a Mars 
lander requires knowledge of the locations of the station and lander, the orientation of Earth and Mars as a function of 
time, the orbits of Earth and Mars including perturbations by the other planets and asteroids, and possible variations 
in the underlying theory of gravity. 

Later simulations were done for a landed transponder on Phobos. The simulations also estimated additional 
body-specific parameters. In particular, the set of Phobos-specific parameters included Phobos' dynamic parameters 
(from integrated Phobos/Deimos ephemeris): Phobos' initial state (position and velocity). Mars' GM (separate from 
planetary list) and tidal lag angle (related to k2/Q), Phobos' GM, J2, and C22, Deimos' initial state (position and 
velocity), and Deimos' GM [72]. We also included Phobos' geometric parameters, namely bias for lander, lander 
coordinates, libration in longitude amplitude (geometric), and libration in latitude amplitude and phase. 

For the simulations using a spacecraft lander on Phobos, the sensitivity to gravitational physics parameters comes 
from the 1-yr orbit of the earth-moon system and the 2-yr-orbit of the Mars system. We did not consider the sensitivity 
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to gravitational physics parameters of the 8 hr Phobos orbit around Mars. Tracking of a lander on Phobos is used as 
a way to follow the motion of more massive Mars through the solar system. 

D. Details of simulation and results 

1. Lander on Mars 

For the Mars lander laser ranging mission, simulated measurements between a Mars lander and an Earth station 
were formed based on a nominal position of the lander and tracking station. Since range measurements within 
a single hour- long observation pass are expected to be correlated (due to pre-pass instrument calibrations), only a 
single "normal" range point was uscid for ciach day. These measurements were then included with the current planetary 
ephemeris data set to estimate parameter uncertainties, with variations depending on assumptions of data accuracy 
and duration. The position of the lander was always estimated from the simulated laser range data along with the 
orbit and Mars rotation parameters. The position of the Earth station was assumed to be known from GPS or another 
surveying technique. 

The orientation of the Earth as a function of time is monitored to the mm level by combinations of VLBI, GPS, and 

laser ranging data. The orientation of Mars as a function of time is modeled in series form with certain parameters 
estimated to fit Mars radio range and Doppler data as part of the planetary ephemeris fit [14, 59]. 

Asteroids have a significant eflfect on the orbit of Mars. For the JPL planetary ephemeris, 67 asteroid mass 
parameters with the largest effect on Mars were individually modeled while 230 other asteroids with smaller eflfects were 
modeled by grouping into three taxonomy classes. In order to investigate the sensitivity of the simulated measurements 
to asteroid modeling, solutions were formed based on allowing various number of asteroid mass parameters to be 
completely unconstrained and requiring estimation from the laser ranging data. 

TABLE III: Estimated uncertainties for parameters of interest as a function of number of estimated asteroid mass parameters, 
for a 36 month Mars lander mission with 1 mm laser ranging once per day with 2° SEP cut-off. 



Rclativistic Effect 
to be studied by PLR 


Number of asteroid GMs 


11 


36 


67 


The Eddington parameter 7 

The Eddington parameter /3 
Strong Equivalence Principle, r/ 
Solar oblatcncss, J 2 
Parameter G/G, yr~^ 


7.8 X 10"'' 

6.9 X lO"-""' 
4.3 X 10"'" 
1.6 X 10"** 
2.6 X lO"^'^ 


1.1 X 10"' 

9.7 X 10"" 

8.8 X 10"" 

2.5 X 10"** 

2.6 X 10"^^ 


1.4 X IQ-' 

1.7 X 10""* 

3.1 X 10""* 

3.2 X 10"* 

2.8 X 10"^'^ 



TABLE IV: Estimated uncertainties for parameters of interest as a function of Mars lander mission duration, with 1 mm laser 
ranging once per day with 2° SEP cut-off and 67 asteroid mass parameters estimated. 



Relativistic Effect 

to be studied by PLR 


Mission duration, months 


18 


36 


72 


The Eddington parameter 7 
The Eddington parameter f5 
Strong Equivalence Principle, ri 
Solar oblateness, J2 
Parameter G/G, yr~^ 


3.1 X 10"' 
4.3 X 10"* 

8.8 X IQ-* 

6.9 X 10"* 
1.7 X 10"" 


1.4 X 10"' 

1.7 X 10"** 

3.1 X IQ-* 

3.2 X 10"* 

2.8 X 10"" 


7.8 X 10"** 
8.6 X 10"^ 
8.5 X 10"^ 
2.1 X 10"* 
1.0 X 10"" 



Table III shows the estimated uncertainties in the PPN 7 and /3 parameters, the second zonal harmonic of the solar 
gravity field, and the rate of change of the gravitational constant G. The estimated uncertainties are clearly dependent 
on the number of asteroid mass parameters that must be estimated. The parameters of most interest, PPN 7 and 
G/G, are less dependent on the asteroid modeling. Since significant improvement in asteroid mass parameters may 
be achieved through such observations and the main parameters of interest for the simulations are not too strongly 
dependent on the asteroid modeling, the estimation of 67 main asteroid mass parameters with the laser ranging data 
is considered to be sufficient for this study. Asteroid mass parameters can also be estimated using observations of the 
deflection of smaller asteroids that happen to pass near larger asteroids. The number of asteroids observed in this 
way is increasing [73]. Also, asteroid masses obtained with other techniques arc a benefit. 

The potential science return, in terms of uncertainties in estimated parameters, is also dependent on mission 
duration. Table IV shows the estimated uncertainties in key parameters as a function of mission duration. This 
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is particularly true for the PPN 7 parameter, which is most strongly measured at solar conjunctions, occurring at 
the Earth-Mars synodic period of 780 days. For a nominal Type-II trajectory,^ landing on Phobos is possible a few 
months before conjunction. A second solar conjunction requires science mission duration of about 3 years, while a 
mission lifetime of 6 years allows observations through three conjunctions. 

2. Lander on Phobos 

Simulations were performed for a lander on Phobos. A Phobos lander avoids issues of dust on the Mars lander. 
For the Phobos lander the orbit and orientation of Phobos must be estimatcid along with the planetary dynamic 
parameters. A numerically integrated orbit of Phobos has been fit to a variety of data by [72]. Sensitivities with 
respect to the initial Phobos orbital elements and Mars tidal parameter were obtained and corrections estimated in 
the overall least-squares fit. To account for possible unmodeled orbital effects, a correction to the orbital elements 
was estimated each year. Phobos rotation was modeled based on a series by [72, 74]. Corrections to the libration 
amplitudes in longitude and in latitude were estimated to account for necessary updates to the model series. Because 
the Phobos orbit and rotation have periods of 8 hours, the relevant parameters are mostly uncoupled with the Mars 
orbital and gravitational parameters of main interest. The simulations performed for Phobos (see Table I) showed 
essentially the same results for gravitational physics parameters as shown in Table IV for the Mars lander case. 

Currently the orbit of Mars is modeled to cm-level accuracy to fit radio range measurements with meter-level 
accuracy, while the orbit of Phobos is modeled to meter-level accuracy to fit spacecraft Dopplcr data for Phobos 
fly-bys from Viking, Mars Express, etc. at the 100 m level. A model of Phobos' rotation, intended to be precise to the 
5 cm level [74], has been calculated analytically. For processing mm-accurate laser ranges, the modeling for Phobos 
will be updated building on lunar experience. The rotation state will be explicitly integrated, and forces affecting 
the orbit of Phobos upgraded to include more accurate tide and solar pressure models. Our PLR simulations include 
estimates of physical parameters needed to achieve the required modeling accuracy, such as the tidal dissipation, 
Phobos J2 and C22, and libration amplitude and phase, and perturbations by asteroids on the Martian orbit. (Note 
that estimated uncertainties for parameters of interest as a function of number of estimated asteroid mass parameters, 
for a 36 month Phobos lander mission with 1 mm laser ranging once per day with 2° SEP cut-off produced results 
nearly identical to those for Mars, shown in Table III.) 

Table I shows the results of parameter estimation arising from simulated laser ranging to Phobos over 1.5-6 years of 
operation based on daily 1 mm accuracy range points. Estimated parameters include orbital elements and masses of 
the planets and Phobos, up to 67 individual asteroids, and densities for 3 classes incorporating another 230 asteroids. 
Based on 3 yr simulations, PLR will be able to measure the PPN parameters 7 and /3 with accuracy levels of 2 x 10~^ 
and 2 x 10^^, respectively, to measure the solar oblateness coefficient J2 with an accuracy of 3 x 10~^, to measure the 
fractional time-rate-of-change of the gravitational constant to 3 x 10~^^ yr"^, and to test the gravitational inverse 
square law to an accuracy of 2 x 10~^^ at -^1.5 AU distances and the SEP parameter to 3 x 10~^. Thus, PLR could 
directly verify important predictions of modern theories with an unprecedented accuracy, enabling strong advances 
in the tests of general relativity. 

An initial set of simulations were performed for a single conjunction as functions of the data accuracy and angular 
separation of Mars from the Sun as viewed from Earth. The primary parameter of interest was the PPN 7, which due 
to the time delay sensitivity, is not as dependent on orbital uncertainties as other parameters of interest. The results 
for the estimated uncertainty in PPN 7 are shown in Figure 5, which are similar to previous results by [23]. Since the 
main science goal is to estimate PPN 7 with an uncertainty of two parts in 10'', these results drive the instrument 
accuracy to 1 mm and require observations within 2° of the Sun (as viewed from Earth). 

3. Lander on Mercury 

Mercury was also investigated as a possible location for the laser ranging lander. With its short orbital period. Mer- 
cury provides many more solar conjunction opportunities than for a comparable mission duration at Mars. Mercury is 
also more sensitive to solar oblateness, and potentially less sensitive to unknown asteroid mass parameters since it is 
farther from the asteroid belt. Simulations were performed along the same lines as the Mars and Phobos simulations. 
Table V shows estimated parameter uncertainties of an 18-month long Mercury surface duration for three different 
sets of estimated asteroid mass parameters. The sensitivity to asteroid modeling is higher than might be expected, 



^ For details, see http://www2.jpl.nasa.gov/basics/bsf4-l.php 
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FIG. 5: Estimated uncertainty in PPN 7 as a function of data accuracy and data cut-off with angular separation from the Sun 
as viewed from Earth. 

probably because Earth is perturbed at a measurable level by the larger asteroids, and the laser ranging is between 
Mercury and Earth. Table VI gives results as a function of mission duration. As expected the estimated uncertainties 
are smaller than in the comparable Mars/Phobos cases, although for the primary parameters of interest, PPN 7 
is improved only by a factor of 3 for the nominal 36-month duration case, and G /G is not significantly improved 
(compare with Table IV). 

TABLE V: Estimated uncertainties for parameters of interest as a function of number of estimated asteroid mass parameters, 
for an 18 month Mercury lander mission with 1 mm laser ranging once per day with 2° SEP cut-off. 



Relativistic Effect 
to be studied by PLR 


Number of asteroid GM s 


11 


36 


67 


The Eddington parameter 7 
The Eddington parameter /3 
Strong Equivalence Principle, r) 
Solar oblateness, J2 
Parameter G/G, yr~^ 


6.7 X lO"** 
4.2 X 10"^ 
2.6 X lO"'^ 
4.4 X 10^^" 
1.2 X 10"" 


8.0 X 10"** 

7.3 X 10"^ 

5.1 X 10"^ 
7.9 X 10"" 

1.4 X 10"" 


9.8 X 10"** 
1.5 X 10"^ 
1.5 X 10"'* 
1.7 X lO"'' 
1.7 X 10"" 



TABLE VI: Estimated uncertainties for parameters of interest as a function of Mercury lander mission duration, with 1 mm 
laser ranging once per day with 2° SEP cut-off and 67 asteroid mass parameters estimated. 



Relativistic Effect 
to be studied by PLR 


Mission duration, months 


18 


36 


72 


The Eddington parameter 7 
The Eddington parameter /3 
Strong Equivalence Principle, 77 
Solar oblateness, J2 
Parameter G/G, yr"^ 


9.8 X 10"** 
1.5 X lO"--' 
1.5 X lO"'' 
1.7 X 10"^ 
1.7 X 10"" 


4.6 X 10"** 
6.6 X 10"'' 
3.4 X 10"^ 
6.9 X 10"" 
2.8 X lO"^'^ 


2.7 X 10"** 
7.5 X 10"^ 
6.5 X 10"*^ 
7.4 X 10"" 
9.1 X lO"^** 



E. Towards PLR Dynamical Model and Secondary Science Objectives 

The development of the PLR data analysis techniques will build on the successful ongoing lunar and planetary data 
analysis. Analysis of PLR data requires modeling the orbital and rotational motions of the lander on Phobos and 
tracking stations on Earth to mm accuracy (see Fig. 6). Models for LLR are currently at cm accuracy for processing 
range measurements to retrorefiectors on the moon. This modeling experience gives confidence that modeling of 
Phobos laser ranging can reach the mm level. The LLR data and models are used to determine the orbit and 
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rotation of the moon, including sensitivity to lunar interior parameters, and establish the state-of-the art for several 
fundamental physics parameters such as the upper bound on the violation of the Strong Equivalence Principle (SEP) 
parameter t], the PPN parameter /3, the value of the temporal variation in the gravitational constant, and the test of 
the inverse-square law of gravity [8] . 




FIG. 6: Precision laser ranging to Phobos could measure the distance between an observatory on the Earth and a terminal on 
the surface of Phobos to accuracy of 1 mm in less then 5 min of integration time. 

A high quality dynamical model of Phobos rotating and orbiting Mars while Mars orbits the Sun is necessary for 
analyzing millimeter quality range data to the PLR spacecraft. Some mission concerns and model items deserve 
special mention and there are also opportunities for science in addition to the gravitational physics objectives. These 
are described below. 



1. Phobos Orbital Science 

Prior to landing on Phobos, the spacecraft would orbit near Phobos for several weeks imaging the surface, measuring 
topography and sampling the gravity field (see discussion in Sec. V A). While the primary purpose of this inspection is 
to provide information to enable a successful landing, there would be a science benefit to mapping surface properties. 
A good determination of the gravity field would aid the development of a physical libration model and recovery of 
moments of inertia. A comparison of topography and gravity field would test whether the density is homogeneous. 
Phobos information has been collected by [75]. 

2. Phobos Physical Librations 

The orientation of Phobos changes with time. A high accuracy physical libration model will be necessary for 
successful range data analysis. The largest departure from uniform rotation has an anomalistic period with an 
amplitude of ~ 1° predicted from the shape. A 1° amplitude would displace a lander at a surface point opposite 
Mars by 230 m. An amplitude of 0.78° ±0.4° has been reported [76]. Measurement of the amplitude will provide the 
moment of inertia combination (B — A)/C, where the moments of inertia satisfy A < B < C . In combination with 
Phobos' gravity coefficient C22 from the mission's orbital inspection, it will be possible to determine the moment of 
inertia C. Another physical libration effect, the equator of Phobos tilts a few seconds of arc with respect to the orbit 
plane. This tilt is sensitive to (C — A)/B\ the combination J2 + 2C22 gives (C — A) while (J2 — 2C22) gives (C — B). 
A theory of Phobos physical librations is given by [74, 77] . Numerically integrated physical librations would need to 
be developed for millimeter data analysis. 

3. Phobos Orbit 

At a first approximation Phobos can be described as orbiting Mars with a precessing slightly elliptical (e — 0.015), 
slightly inclined {i — 1.1°) orbit with semimajor axis a = 9375 km [72]. The orbit period is 7.654 hr and the mean 
motion n is 3.136 revolutions/day, faster than Mars' rotation. The major cause of the precession of perifocus and 
node directions is the oblate figure of Mars. The precession periods are 2.26 yr prograde for perifocus and 2.26 yr 
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retrograde for node. Precession periods shorter than the mission duration benefit separation of parameters during 
solutions. Small additional perturbations arise from Mars' gravity field (other than the oblate figure), the Sun and 
Deimos. An orbit theory has been derived by [78]. 

Perturbations from solar radiation pressure, Mars static and variable gravity field, and tides deserve special discus- 
sion below. We also note that while Phobos' orbit plane precesses along Mars' equator plane, that equator plane is 
precessing and nutating with time. The precession of Mars' equator is of interest because, in combination with Mars' 
J2, it provides C/MB? [14, 57, 58, 79, 80]. 

4.. Solar Radiation Pressure 

Solar radiation pressure acting on Phobos displaces the orbit. The acceleration from solar pressure and thermal 
re-radiation depends on the projected Phobos area divided by its mass.^ This effect has been studied for the Earth's 
moon because it causes a perturbation similar to that from a violation of the equivalence principle [55, 81]. Unlike 
the moon, the projected area of Phobos varies about ±20%. A spherical Phobos would be perturbed by 6.4 mm (vs 
3.65 mm for the moon) in orbit radius and about twice that in orbital longitude. The period is the 0.319 d synodic 
period with respect to the Sun. The real non-spherical Phobos would add harmonics, and the elliptical Mars orbit 
would cause a ±19% modulation of the radial amplitude which would vary from 5 to 8 mm with the 687 d Mars orbit 
period. While the effect of Phobos' shape could get complicated in detail, the relatively small size of the perturbation 
suggests that a smooth shape model would suffice for the acceleration model. 

Mars' obliquity is 25.2°. When the Sun is near either of Mars' equinoxes, there will be eclipses of Phobos by Mars 
and also of Mars by Phobos. This can complicate the solar radiation pressure perturbations of Phobos' orbit by 
allowing a small accumulation of perturbation in orbit longitude, perhaps as much as a few centimeters in several 
months. When the Sun is near either solstice, there are no eclipses. Solar radiation pressure and thermal re-radiation 
could be modeled with four parameters: an acceleration parameter proportional to the area/mass ratio, a thermal lag 
parameter, and two shape parameters. It should be possible to model solar pressure and thermal re-radiation effects 
reasonably well and to solve for the model parameters. 

5. Mars Gravity 

The gravity field of Mars has been mapped with increasing accuracy by successive orbiting spacecraft. The spinning 
planet is oblate and the unnormalized J2 is 1.96 x 10~^. The large J2 causes the Phobos orbit node and perifocus 
directions to precess by 159°/yr, retrograde for the node and prograde for the perifocus. While Phobos' eccentricity 
and inclination are small, these precession rates would be determined to very high accuracy by Phobos laser ranging. 
Since the sensitivity to higher-degree harmonics falls off as powers of 1/2.76 (the ratio of the radius of Mars R to 
Phobos' semimajor axis a), spacecraft-derived static fields for Mars should be adequate in higher degrees, but it may 
be necessary to solve for some lower degree gravity coefficients. 

Of particular interest are time variations in Mars' gravity field. Seasonal changes of the polar caps, with exchange 
of CO2 between the poles and variable atmospheric pressure, cause observed variations in the gravity field [14, 79, 82]. 
Seasonal changes in J3 are clearly seen from the gravity variations while possible variation in J2 is closer to the noise. 
Variation in J2 is proportional to variation in the polar moment of inertia, which is independently observed through 
the rotation of Mars determined from radio tracking of landed spacecraft [57, 59] and orbiting spacecraft [14, 80]. J4 
and J5 variations are also reported by [83] . Phobos has a low inclination orbit so it is most affected by the even zonal 
coefficients. A 1 x 10~^ variation of the unnormalized J2 with 687 d period would cause a term with 8 m amplitude 
in Phobos' orbital longitude. The peak-to-peak seasonal variation in J2 predicted by [83] is 7 x 10~^. Using this 
simulation and the rotation information for J2 variations, orbit perturbations with periods of 687 d and 343 d should 
be prominent (one to several meters). Higher harmonics should be present with smaller amplitudes. Seasonal J2 
changes, combined with other even seasonal zonal variations, should be strongly detected by PLR. 

Accurately tracking a lander on Phobos offers great promise for determining effects from the seasonal waxing and 
waning of the Martian polar caps. At some level, these seasonal changes will not exactly repeat from year to year and 
that will lead to unpredictable dynamical perturbations of Phobos' orbit. Seasonal changes and small nonperiodic 
variations in J2 will have to be part of the solution. We do not expect that a small long time scale wander of the 



The values for the area-to-mass ratio for the moon, Phobos and Mars were evaluated to be (j4/M)moon — 1.29 x 10 m^kg ^ for the 
moon, (A/M)phobos ^ 3.3 x 10"* m^kg-^ for Phobos, and (A/M)Mars ^ 5.65 x 10"" m^kg-^ for Mars. 
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orbital longitude of Phobos, which would affect the range observable at short period, would limit the science goals 
of long period in the center of mass motion of Mars and its satellites around the Sun since that center of mass can 
be accurately determined over short time spans. Simulations, with the Phobos orbit independently determined every 
year, successfully recovered the gravitational physics parameters. 

6. Tides on Mars 

Solid-body tides are raised on Mars by both Phobos and the; Sun. Tides from Phobos are about 1 mm in height 
while Sun-raised tides are ~ 1 cm. Gravitational variations from tides on Mars can perturb the orbit of Phobos. 
Sun-raised zonal tides have long periods, 687 d and 343 d, from orbit eccentricity and obliquity respectively. The 
corresponding Phobos orbit perturbations are about 1 m in longitude. There are also smaller tidal signatures at 1/3 
and other fractions of the 687 d Martian year. The periods are the same as the polar cap J2 variations, but the phases 
are known for tides. A tide model should be included for accurate dynamics. 

7. Phobos Tidal Acceleration 

Phobos shows a large acceleration in orbital longitude. Recent fits by Bills et al. [84], Lainey et al. [85], and 
Jacobson [72] give accelerations in the forward orbital longitude adn/dt of ~ 416 m/yr^. This secular acceleration 
would easily be detected by PLR giving refined accuracy. The apparent cause of this acceleration is Phobos-raised 
tides on Mars perturbing Phobos. The tidal bulge raised by Phobos is behind (in time and longitude) Phobos' position 
sapping energy from the orbit, which consequently shrinks by 3.8 cm/yr. Phobos will eventually impact Mars [86]. 
The most important of the tidal components for the secular acceleration should be the second-degree M2 tide of 
period 5.55 hr on Mars. The small eccentricity (0.015) and inclination (1.1°) tend to reduce the influence of other 
degree 2 tides by ^ 3 orders of magnitude or more. The influence of tides of higher degree fall off as even powers 
of R/a = 1/2.76, about an order of magnitude per degree. So the degree 3 tide of 3.7 hr period on Mars is a small 
contributor to the tidal secular acceleration. 

What do tides on Phobos do? Based on evolutionary arguments, Yoder [87] has placed an upper limit on Phobos' 
k2/Q of 2 X 10~^, which would make dissipation in Phobos a minor contributor, of order 10""^ relative to the overall 
tidal acceleration of Phobos. Periodic tidal displacements on Phobos might reach 1 mm. 

8. Impacts 

Mars and other solar system bodies are subject to bombardment by meteoroids. Small new craters are observed 
on Mars [88]. The probability of impact on a body depends on the body's projected area while the recoil from the 
impulse is inversely proportional to the body's mass. The area to mass ratio for Phobos is ~600 times larger than 
that for Mars, so impacts are considered here. A 1 kg projectile impacting Phobos at 10 km/s can cause a mean 

motion change corresponding to a An up to 0.1 mm/yr along the orbit. Inner solar system impact rates would make a 
1 kg or larger impact on Phobos a rare event occurring less than once per century [89-92] . Impacts are not a concern 
for the dynamics of the PLR mission. 

9. Asteroid Masses 

Mars and Earth are perturbed by minor planets in the asteroid belt. Mars is adjacent to the belt and is much 
more strongly influenced than the Earth. Both analytical and numerical techniques have shown the importance of 
asteroid perturbations [14, 42, 93, 94] . The three largest asteroids, 1 Ceres, 2 Pallas and 4 Vesta, are each capable of 
perturbing Mars by more than one kilometer on time scales of several decades through long-period near resonances. 
Shorter period perturbations with smaller amplitudes are also effective. For periods of a decade or less, amplitudes 
are > 100 m from Ceres and Vesta, somewhat smaller from Pallas, and at least 10 m from four other bodies. There 
are estimated to be ~ 67 bodies which can perturb Mars with amplitudes > 1 m for periods of 10 yr or less and the 
number is hundreds for several decades. The total perturbation is larger than any single amplitude. Knopoliv et al. 
found 67 bodies which perturb the Earth-Mars range by > 100 m long term or > 4 m short term. 

At a millimeter level Mars is affected by a very large number of small bodies, too many to reliably separate all 
of the masses during a solution. The simulations of this paper solved for 67 asteroid masses. The asteroid masses 
affect the simulated solutions for several parameters of interest including the equivalence principle. The acceleration 
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in orbital longitude from solar GM change is difBcult to distinguish from periodicities much longer than the data 
span. Existing Earth-Mars radio range data can be fit well by solving for 20 or more asteroid masses [14, 42, 60, 80], 
but the asteroids cause the orbit error to grow rapidly for extrapolations into the future. 

Apart from improving the dynamics of the solar system, asteroid masses are of scientific interest. Asteroids have 
different compositions with physical properties grouped into multiple taxonomic classes indicating a diversity of origin 
and evolution histories. Their densities are part of the physical properties picture. In addition to fits of radio ranges, 
existing determinations of asteroid masses come from asteroid-asteroid perturbations observed by optical astrometry 
[73], observations of binary asteroids, several spacecraft flybys and one orbiter. Masses derived from very accurate 
laser ranges to Mars or Phobos would be sensitive to the more massive asteroids and those that can pass through the 
inner asteroid belt. The future Gaia astrometric mission is expected to determine masses from perturbations of one 
asteroid passing another [95]. 

V. TECHNICAL OVERVIEW OF PLR MISSION, EXPERIMENT, AND METHODOLOGY 

A. Mission Description 

The Phobos Laser Ranging mission is based on deploying a laser ranging transponder instrument capable of sup- 
porting measurements of distances to Earth with 1 mm accuracy given daily hour- long tracking passes [24]. In our 
studies, laser ranging in the optical regime was selected rather than radio ranging. While theoretically it may be 
possible to achieve mm class accuracy using a dual Ka and X band arrangement, laser ranging is much closer to 
the state of the art. Furthermore, ground-based radio antenna calibration would be a serious concern-currently the 
DSN does not calibrate to even meter accuracy in range. A pulsed, photon counting laser system was selected rather 
than a coherent system. The number of received signal photons necessary for coherent detection was not considered 
attainable at standard laser power levels at interplanetary ranges. 




FIG. 7: (a) PLR spacecraft in deployed configuration, (b) Phobos landing trajectory concept. 

The PLR spacecraft configuration is fairly simple. The box-like main structure would have four legs rigidly attached 
for landing on Phobos (Fig. 7). Given the anticipated characteristics of the spacecraft, it will easily fit within the 
fairing of an Atlas V-401 launch vehicle, the smallest Mars-capable vehicle under NASA contract for 2015 and beyond. 
After launch the four small solar panels, high-gain antenna, and laser ranging instrument will be deployed. 

Following a 7-month cruise, when the spacecraft arrives at Mars, a propulsive maneuver will place it into a highly 
elliptical orbit about Mars. This will be followed by a maneuver at apoapsis to raise the periapsis and match the 
plane of Phobos' orbit, and then a maneuver at peripasis to match the Phobos orbit period (Fig. 7b). Phobos will 
be observed from orbit for one month with a simple camera to improve the Phobos ephemeris and identify a specific 
landing location and map gravity field. (The PLR ranging system could be designed to operate in altimeter mode 
before landing, which may help to improve the determination of the Phobos' gravity field.) Landing on Phobos will 
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be done semi-autonomously using a laser altimeter (could be separate from PLR transponder) and a feature detection 
camera similar to those used on the NEAR spacecraft/ 
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FIG. 8: (a) Sequence of laser pulse time-tags used to derive Earth-Phobos round-trip light-time, (b) A block diagram of the 
laser ranging instrument and Earth instrument for timing transmitted and received pulses on a common detector (see more 
details on PLR instrument in Fig. 10). 



After landing, the spacecraft will have approximately three sun and Earth view periods of > 3 hours, giving more 
than 5 hours with elevations > 45° every Earth day. One of these view periods each day will be used for the hour-long 
laser ranging experiments. Once each week, a second sunlit period will be used to send the instrument data to Earth 
via the high-gain antenna. Nominal science operations would continue for 3 years, covering two solar conjunctions. 
Since there are no expendables (e.g. propellant) during science operations, further observations could occur in an 
extended mission phase. 



B. Measurement Description 



The gravity parameter estimation is based on monitoring the light-time between Earth and Phobos. Unlike lunar 
laser ranging which measures the round-trip light-time from Earth to the moon and back, PLR will use a combination 
of one-way light-times with asynchronous operation of the laser [13, 20]. The advantages of the asynchronous scheme 
include (for details on laser transponder operation see [13]): (1) noise immunity: the remote station does not transmit 
pulses in response to numerous background photon detections; (2) stable laser performance: steady pulse generation 
results in thermally stable conditions; (3) the measurement rates for an "asynchronous transponder" are substantially 
higher compared to an "echo" one, if the probability of detection at one or both ends of the link are small; and (4) 
multiple Earth stations [10] can range to the remote unit simultaneously without placing conflicting demands on the 
remote transmitter. The Earth stations will measure the times at which laser pulses are directed toward Phobos and 
the times of received photons. Similarly, the laser ranging instrument on Phobos will measure the times of pulses 
transmitted towards Earth and the times of received photons from Earth. Scattered light will be rejected to first 
order by using a known laser pulse rate and filtering the measured times to find those separated by the signal rate. 
Timing will be done to an accuracy of 50 ps per photon (though not explicitly synchronized), averaging down to < 3 
ps for a typical observation. For the Earth station, the timing accuracy must be kept over the round-trip light time 
of up to ^3000 s. The Phobos instrument accuracy need only be good over the few minutes it takes to collect enough 
photons to produce one-millimeter statistics. The timing scheme is indicated in Fig. 8a. 

Both the Earth observatory and laser ranging instrument include a telescope and optics to have incoming light 
and a small fraction of outgoing light illuminate a photon sensitive detector array that records single photon arrivals 
(Fig. 8b). The power received by the telescope depends directly on the telescope's collecting area and inversely on 



See more details on NEAR spacecrfat at http://www.msss.eom/small_bodies/near_iiew/index.html#top 
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the returning spot area. At the detector, both a pinhole restricting the field of view and a narrow-band-pass filter 
will be included to reduce background light. 

If the combined timing uncertainty per photon is 50 ps (from laser pulse width, detector response, timing electronics, 
clock), then one needs (50/3.3)^ = 225 photon events in order to establish the one-way distance to 1 mm (3.3 ps) 
precision. At the estimated rates in the most demanding scenario, this is achieved in a few minutes of ranging. 
However, it is assumed that all other parameters (beam pointing, divergence) are satisfied, which will not be true all 
of the time, so hour-long tracking passes are used to ensure that enough valid measurements are obtained. 

TABLE VII: Signal and noise estimated for laser ranging from Earth to Phobos and back-most demanding scenario. Note: "tx" 
and "rx" stand for transmit and received signal modes. Also, the following notations are used pW= 10~^^W, aW= 10~'^*W, 
pj= 10~'2j, aJ= 10-'*J. 



Link parameters 


Earth to Phobos 


Phobos to Earth 


Input parameters: 


wavelength 


532 nm 


1064 nm 


tx power 


3 W 


0.25 W 


tx efficiency 


0.5 


0.5 


tx beam divergence 


25 fjxeud 


160 fjxad. 


tx pointing loss 


-2 dB 


-2 dB 


tx atmospheric loss 


-3 dB 


dB 


tx pulse frequency 


1 kHz 


1 kHz 


rx atmospheric loss 


dB 


-4.3 dB 


rx diameter 


0.1 m 


1 m 


rx efficiency 


0.3 


0.3 


rx field of view 


240 /xrad 


20 /xrad 


rx detector efficiency 


0.4 


0.4 


background 


65 W/m'^/sr//im 


20 W/m^/sr/^im 


scattered light irradiance 


130 W/mVsr/^m 


130 W/mVsr/^im 


Earth's sky irradiance 


W/m7sr//jm 


100 W/m7sr//im 


bandpass FWHM 


0.2 nm 


0.2 nm 


range 


2.6 AU 


2.6 AU 


Derived parameters: 


photon energy 


0.37 aJ 


0.18 aJ 


space loss 


-379.3 dB 


-373.2 dB 


rx signal power 


36.1 aW 


37.8 aW 


planet background power 


0.16 pW 


1.2 pW 


scattered light power 


17 pW 


2.5 W 


sky radiance power 


pW 


6.3 pW 


Summary results: 


incident signal power 


10.8 aW 


1.13 aW 


incident noise power 


5.1 pW 


3.0 pW 


SNR 


-56.8 dB 


-64.2 dB 


detected signal rate 


11.6 Hz 


2.4 Hz 


detected noise rate 


5.5 MHz 


6.5 MHz 


10 ns window data volume 


190 Mb/day 


222 Mb/day 



The laser ranging instrument on Phobos will transmit pulses at a wavelength of 1064-nm for greater power efficiency 
and lower atmospheric scattered light at the Earth receiver during near sun pointing operations. A laser pulse rate 
of 1 kHz will be used with 12 ps pulse-widths and power of 0.25 W. We selected a divergence of 160 /irad to cover 
the entire Earth at distances greater than 1 AU. At closer ranges, the instrument must point to a specific area on the 
Earth. For aim aperture on the Earth station at maximum Earth-Mars range, the average detection rate will be 
1.2 photons per second. 

Transmission from the Earth station will be at a wavelength of 532 nm and detected by the instrument on Phobos 
with an efficient, low noise silicon-based single photon detector. The beam transmitted from Earth will have an uplink 
divergence of 25 /^rad to cover potential pointing and atmospheric seeing variations. Assuming a laser power of 3 W 
(3 mJ/pulse at 1 kHz repetition rate) and a 12 cm aperture on Phobos, the uplink will deliver 3 detected-photon/s 
at the maximum Earth-Mars range. 

Background light on Earth represents a serious challenge for a detector system capable of single-photon signal 
detection. Every photon is identical as seen by the detector, so that signal photons may only be differentiated 
from background photons via time-domain characteristics. Each square arcsec of the illuminated Earth delivers 2000 
photons per second into a 0.01 m^ aperture employing a 0.2 nm bandpass filter and 10% overall detection efficiency. 
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Assuming the detector has a spatial resohition of 1 arcsec, this means that background hght will typically dominate 
the signal in photon count. It becomes clear why echo-transponding is less favorable in these conditions. However, the 
laser pulses — even if only occasionally detected — will form a regular temporal pattern unlike the background photons, 
and will therefore be separable [13]. 

Refractive delay of the earth atmosphere imposes a '^2 meter propagation delay, which is routinely removed at the 
few mm level via surface measurements of pressure, temperature, and humidity [96, 97]. 

The background event rates-due to scattered sunlight, illuminated body in field-of-view, and the radiance of the 
Earth atmosphere-will be much greater than signal event rates. At Phobos the laser ranging instrument will search 
for pulses in a 1 kHz repeat pattern following the tracking passes, and select only pulses within a narrow window 
about the times of the expected pulses from Earth for transmission. At Earth the photon timings will be matched 
against predicted pulse rate and Phobos-Earth range to select correct pulses. Adjusting the pulse repetition pattern 
on Earth to be precisely 1 kHz at Phobos will facilitate this process. Table VII summarizes a conservative estimate 
of signal and noise levels at the most distant separation — at conjunction when the sky is also at its brightest (note, 
"tx" and "rx" in the Table VII stand for transmit and received). 



C. System Description 

1. Ground Segment System 

One or more Earth observatories will be equipped similar to current lunar laser ranging stations [10] and be outfitted 
with membrane-type sunlight rejection filters for accurate measurements when Mars is near solar conjunction. The 
ground stations will use lasers very similar to those on Phobos, with the exception that the frequency will be doubled 
and the pulse repetition rate will be continuously tunable so that the pulse arrival rate at Phobos is exactly 1 kHz. 
The detector used in the ground stations will be an InGaAsP hybrid photodiode with a thinned photocathode to 
reduce timing jitter to less than 40 ps. A solar rejection filter with passband at 532 and 1064 nm must cover the 
entire ~1 m aperture of the telescope, which is a demanding task, but one that has been previously demonstrated. 




2. Payload 

The instrument design is evolved from optical communications 
terminals under development for interplanetary missions. Proto- 
types of the communications terminals have been tested in air- 
planes. For the Phobos laser ranging science objectives, lower pulse 
rates are acceptable but higher timing resolution is required. The 
PLR instrument design is based on existing laboratory lasers and 
array photo-detectors which will be ruggedized to operate in space 
[25]. The instrument is shown in Fig. 9 and a functional block 
diagram shown in Fig. 10. The instrument component mass and 
power are listed in Table VIII. 

The instrument consists of a gimbaled 12 cm telescope, a 250 
mW Nd:YAG laser operating at 1064 nm, a silicon Geiger-mode 

(GM) avalanche photodiode detector (APD) array in 8x8 format, ^^^^^ 1^^^^ ^^^^.^^ instrument 

a multiplexed ASIC timing system to go with the detector, a clock, telescope. 

and an acquisition/pointing CCD camera with thermo-electric cooler. A small actuated mirror will impart a point- 
ahead angle to the beam with respect to the receiver direction of up to 330 /irad, which will be monitored via the 
centroid of a sample of the transmit beam imaged onto the CCD. A small fraction of the outgoing laser light will be 
reflected back from a reference mirror to the single-photon detector array for timing transmitted pulses. The laser 
will derive from a Semiconductor Saturable Absorber Mirror (SESAM) fiber seed laser, delivering 12 ps pulses at a 
50 MHz repetition rate, followed by a pulse-picker that will accept only every 50,000th pulse to form a 1 kHz train. 
A laser diode-pumped regenerative amplifier will then deposit ^250 fiJ of energy into each pulse for transmission to 
Earth. 

The instrument electronics includes timing, pointing control, data storage, and a data post-processor. After ac- 
quiring photon times for the one hour tracking pass, the post-processor will look for the 1 kHz transmit pulse pattern 
from Earth, delete times outside a 10 ns window about that estimated rate, and format times into differences for 
reducing data volume. Edited times will be collected for transmission to Earth via radio, with a data volume of 5 Mb 
per day. 
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FIG. 10: PLR instrument block diagram. 



3. Flight System 

The PLR spacecraft includes all equipment necessary to deliver and support the laser ranging instrument on 
Phobos. The spacecraft structure is a simple box-like construction to accommodate the payload, spacecraft electronics, 
propulsion, attitude control, telecommunications, thermal, and power subsystems. The mass by subsystem is given 
by Table IX. There are several operating modes; a typical power by subsystem is also given in Table IX. With the 
Atlas V 401 as the launch vehicle, mass and volume capabilities are much greater than necessary, so subsystems have 
been optimized to minimize cost rather than mass. Selected redundancy has been included in subsystem design to 
reduce risk where cost effective. 

Attitude control will be done using thrusters since there is no stringent pointing requirement for the spacecraft. 
Attitude control sensors include standard star trackers, inertial measurement units, coarse sun sensors, optical nav- 
igation camera, terminal descent camera, terminal descent laser altimeter (i.e., a NEAR-type altimeter mentioned 
earlier), and electronics for driving the laser range instrument pointing gimbal. 

The command and data handling subsystem is based on the architecture being developed for the Mars Science 
Laboratory (MSL). It includes RAD750 processors, interfaces to the attitude control, telecommunication, power, and 
payload subsystems, with sufhcient memory storage to accumulate up to two weeks of instrument data. 

The power subsystem includes 5.3 m^ solar panels and batteries to operate in absence of sun-light. The battery 
size is dictated by the propulsion system operation during Mars orbit insertion. 

Propulsion will be done using a mono-propellant to minimize cost of tanks and thrusters. Standard thrusters of 
300 N(4), 90 N(4), 4.5 N(4), and 0.9 N(12) will be used for Mars orbiter insertion, attitude control, and landing on 
Phobos. Four propellant tanks will contain adequate propellant with margin. After landing the remaining propellant 
will be vented. 

Thermal control will be done using standard multi-layer insulation and heaters. Spacecraft and electronics will be 
located in a thermally isolated box that will be kept warm during dark periods on Phobos. The propellant tanks will 
be kept warm until landing and uncontrolled afterwards. 

Telecommunications will utilize standard X-band transponders and traveling-wave tube amplifiers, with a deployed 
pointable 1.5-m diameter high-gain antenna for sending instrument data to Earth, and two low-gain antennas used 
for contacts during cruise. 
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TABLE VIII: Laser ranging instrument mass/power. 



Component | Mass, kg | Power, W 
Subassembly: Optics 



Primary mirror 


0.665 




Secondary mirror 


0.018 




Sun filter 


0.475 




Mirror doublets 


0.009 




Fold mirrors 


0.008 




Beam splitters 


0.062 




Subassembly: Electronics 


Power conditioning 


3.000 


5.25 


Controller (FPGA) 


3.000 


7 


Seed laser 


0.500 


0.25 


Pump laser 


6.000 


10.5 


Timing storage 


0.500 


2.2 


GM-APD readout 


1.000 


3 


CCD readout 


1.000 


3 


Gimbal drive 


1.000 


10 


Point-ahead drive 


0.250 


1 


Shutter drive 


0.500 


3 


Clock (USO) 


1.000 


2 


Thermal control 


0.500 


2 


Subassembly: Mechanical 


Telescope gimbal 


4.500 




Point-ahead actuator 


0.065 




Shutter activator 


0.200 




Primary structure 


4.800 




Secondary structure 


0.580 




Subassembly: Thermal 


Heat pipe 


4.000 


1 


Insulation 


0.500 




Survival heaters 


0.800 




Thermal switches 


0.400 




Cooler for CCD 


0.300 


1 


Thermometers 


1.000 




Thermal strap 


2.400 




Fine control heaters 


0.040 


1 


Fine control thermometers 


0.200 




Total: 


39.272 


52.200 



TABLE IX: Flight System Mass and Power Estimates. 



Subsystem 


Mass, kg 


Peak Power, W 


Peak Power Mode 


Laser Ranging Instrument 


39.3 


50 


Science 


Attitude Control 


29.4 


48 


Landing 


Command and Data 


22.0 


47 


Science 


Power 


71.4 


71 


Orbit Insertion 


Propulsion 


170.6 


217 


Orbit Insertion 


Structure/Mechanisms 


203.6 







Launch Adaptor 


14.8 







Cabling 


40.4 







Telecom 


21.0 


64 


Data Ti'aiismit 


Tlienuai 


38. 1 


n>r) 


(,'rui.so 


Spacecraft Total (dry) 


650.9 


619 


Orbit Insertion 


Contingency 


280.1 


266 




Spacecraft -|- Contingency 


931.0 


885 




Propellant 


1216.2 






Spacecraft Total (wet): 


2147.2 
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4- Mission Operations System 

Science operations will be done for one-hour tracking passes each Earth day. The spacecraft will have three view 

periods each day to select for operations. The tracking pass time will be selected to accommodate a view of the selected 
ground observatory for that day. Observing times will be optimized to have the observatory view Phobos at a high 
elevation, and sample the Phobos orbit by varying the time to alternate between before and after maximum elevation 
at the spacecraft. Times will be selected one week in advance. Once per week the spacecraft will communicate with 
the Deep Space Network to download the previous week of science data and to upload the schedule of tracking times 
for the next week. 



D. Science Data Collection, Analysis, and Archive 

The primary data product from PLR will be a series of timestamps for photon events, each also indicating pixel of 
origin on the detector. The signal-finding algorithms will permit sorting the local fiducial return (launch times) from 
one or more receive pulse trains from stations on the earth. Compressed statistics on background rate for each pixel 
will also be available. To the extent permitted by communications capabilities, full-rate data (including background 
events) will sometimes be included so that it will be possible to verify on-board compression efficacy. Ground stations 
will likewise present timestamps for their photon launch and receive events. Auxiliary data such as meteorological 
conditions at each site, system health parameters, pointing information, etc. will join the total dataset, but will 
constitute a small fraction of the data volume. 

Data processing occurs in two stages. First, the event times are combined into meaningful "normal points," 
comprised of a representative midpoint time and round-trip time constructed from the two one-way measurements. 
During this process, the Phobos clock is calibrated (frequency ofi^set measured), and various detector and timing system 
performances are verified. Second, having a set of normal points from a variety of earth stations, the measurements 
are compared to a sophisticated, parameterized physical model of the solar system, including gravitational dynamics 
and planetary physics influences. The difference between the model and actual ranges is minimized with respect to 
the parameters in the model-many of which are simply initial conditions for solar system bodies. The model for 
gravity allows exploration of departures from general relativity, which is the primary focus of the PLR project. 

VI. CONCLUSIONS 

Interplanetary laser ranging may offer very significant improvements in many areas including fundamental physics, 

planetary research and deep-space navigation. What is critical for the purposes of fundamental physics is that, while 
in free space, ILR allows for a very precise trajectory estimation to an accuracy of less than 1 cm at distances of 
^ 2 AU. With the recent successful laser transponder experiments conducted with MLA and MOLA instruments 
[15-17, 19], ILR rapidly becoming a mature technology. A mm-level ranging precision over interplanetary distances 
is within reach, thus opening a way for significant advances in the tests of gravity on solar system scales [18]. 

Laser ranging to Phobos with mm-level range sensitivity would allow major advances in the tests of fundamental 
gravity [3]. A relativistic time delay of the laser signal passing through the Sun's gravity field will deliver a high 
accuracy determination of the Eddington parameter 7. It is important to improve the accuracy of tests of the 
parameter 7 since any deviation from general relativity would signal new physics. Although, no break-down of the 
PPN formalism is expected, the unprecedented precision of PLR could point out where the effects oc on light 
propagation kick in. We emphasize that PLR offers the opportunity to significantly improve the measurement of the 
time-rate-of-change in the gravitational constant G, and to search for a new long-range interaction via tests of the 
Yukawa forces on interplanetary scales. In addition, PLR will test the EP looking for Jupiter-induced deviations in 
the Earth-Phobos range as these bodies orbit the Sun. 

In addition to the primary astrophysics science, a variety of other scientific advances would be made by the PLR 
mission. For Phobos, these include detailed gravity, topographical, and thermographic maps, unprecedented orbital 
determination acciiracy and knowledge of physical libration characteristics, an upper limit on tidal behavior, and local 
information about the landing site, such as rock distribution, regolith characterization, etc. For Mars, its overall body 
distortion and dissipation due to tides could be measured. 

The future deployment of laser transponders on interplanetary missions will provide new opportunities for highly 
improved tests of gravitational physics. With their anticipated capabilities, interplanetary transponders will lead 
to very robust advances in the tests of fundamental physics and could discover a violation or extension of general 
relativity, or reveal the presence of an additional long range interaction in physical laws. As such, these devices should 
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be used in planning both the next steps in lunar exploration and future interplanetary missions to explore the solar 
system. 
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